Genome composition and GC content influence loci distribution in reduced representation genomic studies.

BACKGROUND: Genomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored. RESULTS: Here, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected. CONCLUSIONS: Our results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.

Mutations and intron polymorphisms in voltage-gated sodium channel genes of different geographic populations of Culex pipiens pallens/Culex pipiens quinquefasciatus in China.

BACKGROUND: Culex pipiens pallens and Culex pipiens quinquefasciatus are the dominant species of Culex mosquitoes in China and important disease vectors. Long-term use of insecticides can cause mutations in the voltage-gated sodium channel (vgsc) gene of mosquitoes, but little is known about the current status and evolutionary origins of vgsc gene in different geographic populations. Therefore, this study aimed to determine the current status of vgsc genes in Cx. p. pallens and Cx. p. quinquefasciatus in China and to investigate the evolutionary inheritance of neighboring downstream introns of the vgsc gene to determine the impact of insecticides on long-term evolution. METHODS: Sampling was conducted from July to September 2021 in representative habitats of 22 provincial-level administrative divisions in China. Genomic DNA was extracted from 1308 mosquitoes, the IIS6 fragment of the vgsc gene on the nerve cell membrane was amplified using polymerase chain reaction, and the sequence was used to evaluate allele frequency and knockdown resistance (kdr) frequency. MEGA 11 was used to construct neighbor-joining (NJ) tree. PopART was used to build a TCS network. RESULTS: There were 6 alleles and 6 genotypes at the L1014 locus, which included the wild-type alleles TTA/L and CTA/L and the mutant alleles TTT/F, TTC/F, TCT/S and TCA/S. The geographic populations with a kdr frequency less than 20.00% were mainly concentrated in the regions north of 38° N, and the geographic populations with a kdr frequency greater than 80.00% were concentrated in the regions south of 30° N. kdr frequency increased with decreasing latitude. And within the same latitude, the frequency of kdr in large cities is relatively high. Mutations were correlated with the number of introns. The mutant allele TCA/S has only one intron, the mutant allele TTT/F has three introns, and the wild-type allele TTA/L has 17 introns. CONCLUSIONS: Cx. p. pallens and Cx. p. quinquefasciatus have developed resistance to insecticides in most regions of China. The neighboring downstream introns of the vgsc gene gradually decreased to one intron with the mutation of the vgsc gene. Mutations may originate from multiple mutation events rather than from a single origin, and populations lacking mutations may be genetically isolated.

Systematic identification and characterization of exon-intron circRNAs.

Exon-intron circRNAs (EIciRNAs) are a circRNA subclass with retained introns. Global features of EIciRNAs remain largely unexplored, mainly owing to the lack of bioinformatic tools. The regulation of intron retention (IR) in EIciRNAs and the associated functionality also require further investigation. We developed a framework, FEICP, which efficiently detected EIciRNAs from high-throughput sequencing (HTS) data. EIciRNAs are distinct from exonic circRNAs (EcircRNAs) in aspects such as with larger length, localization in the nucleus, high tissue specificity, and enrichment mostly in the brain. Deep learning analyses revealed that compared with regular introns, the retained introns of circRNAs (CIRs) are shorter in length, have weaker splice site strength, and have higher GC content. Compared with retained introns in linear RNAs (LIRs), CIRs are more likely to form secondary structures and show greater sequence conservation. CIRs are closer to the 5'-end, whereas LIRs are closer to the 3'-end of transcripts. EIciRNA-generating genes are more actively transcribed and associated with epigenetic marks of gene activation. Computational analyses and genome-wide CRISPR screening revealed that SRSF1 binds to CIRs and inhibits the biogenesis of most EIciRNAs. SRSF1 regulates the biogenesis of EIciLIMK1, which enhances the expression of LIMK1 in cis to boost neuronal differentiation, exemplifying EIciRNA physiological function. Overall, our study has developed the FEICP pipeline to identify EIciRNAs from HTS data, and reveals multiple features of CIRs and EIciRNAs. SRSF1 has been identified to regulate EIciRNA biogenesis. EIciRNAs and the regulation of EIciRNA biogenesis play critical roles in neuronal differentiation.

The regulatory roles of small nucleolar RNAs within their host locus.

Small nucleolar RNAs (snoRNAs) are a class of conserved noncoding RNAs forming complexes with proteins to catalyse site-specific modifications on ribosomal RNA. Besides this canonical role, several snoRNAs are now known to regulate diverse levels of gene expression. While these functions are carried out in trans by mature snoRNAs, evidence has also been emerging of regulatory roles of snoRNAs in cis, either within their genomic locus or as longer transcription intermediates during their maturation. Herein, we review recent findings that snoRNAs can interact in cis with their intron to regulate the expression of their host gene. We also explore the ever-growing diversity of longer host-derived snoRNA extensions and their functional impact across the transcriptome. Finally, we discuss the role of snoRNA duplications into forging these new layers of snoRNA-mediated regulation, as well as their involvement in the genomic imprinting of their host locus.

Intron detention tightly regulates the stemness/differentiation switch in the adult neurogenic niche.

The adult mammalian brain retains some capacity to replenish neurons and glia, holding promise for brain regeneration. Thus, understanding the mechanisms controlling adult neural stem cell (NSC) differentiation is crucial. Paradoxically, adult NSCs in the subependymal zone transcribe genes associated with both multipotency maintenance and neural differentiation, but the mechanism that prevents conflicts in fate decisions due to these opposing transcriptional programmes is unknown. Here we describe intron detention as such control mechanism. In NSCs, while multiple mRNAs from stemness genes are spliced and exported to the cytoplasm, transcripts from differentiation genes remain unspliced and detained in the nucleus, and the opposite is true under neural differentiation conditions. We also show that m6A methylation is the mechanism that releases intron detention and triggers nuclear export, enabling rapid and synchronized responses. m6A RNA methylation operates as an on/off switch for transcripts with antagonistic functions, tightly controlling the timing of NSCs commitment to differentiation.

Post-transcriptional splicing can occur in a slow-moving zone around the gene.

Splicing is the stepwise molecular process by which introns are removed from pre-mRNA and exons are joined together to form mature mRNA sequences. The ordering and spatial distribution of these steps remain controversial, with opposing models suggesting splicing occurs either during or after transcription. We used single-molecule RNA FISH, expansion microscopy, and live-cell imaging to reveal the spatiotemporal distribution of nascent transcripts in mammalian cells. At super-resolution levels, we found that pre-mRNA formed clouds around the transcription site. These clouds indicate the existence of a transcription-site-proximal zone through which RNA move more slowly than in the nucleoplasm. Full-length pre-mRNA undergo continuous splicing as they move through this zone following transcription, suggesting a model in which splicing can occur post-transcriptionally but still within the proximity of the transcription site, thus seeming co-transcriptional by most assays. These results may unify conflicting reports of co-transcriptional versus post-transcriptional splicing.

Multiplexed gene editing in citrus by using a multi-intron containing Cas9 gene.

Several expression systems have been developed in clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein 9 (CRISPR/Cas9) framework allowing for gene editing of disease-associated genes across diverse citrus varieties. In this study, we present a new approach employing a multi-intron containing Cas9 gene plus multiple gRNAs separated with tRNA sequences to target the phytoene desaturase gene in both 'Carrizo' citrange and 'Duncan' grapefruit. Notably, using this unified vector significantly boosted editing efficiency in both citrus varieties, showcasing mutations in all three designated targets. The implementation of this multiplex gene editing system with a multi-intron-containing Cas9 plus a gRNA-tRNA array demonstrates a promising avenue for efficient citrus genome editing, equipping us with potent tools in the ongoing battle against several diseases such as canker and huanglongbing.

A concerted increase in readthrough and intron retention drives transposon expression during aging and senescence.

Aging and senescence are characterized by pervasive transcriptional dysfunction, including increased expression of transposons and introns. Our aim was to elucidate mechanisms behind this increased expression. Most transposons are found within genes and introns, with a large minority being close to genes. This raises the possibility that transcriptional readthrough and intron retention are responsible for age-related changes in transposon expression rather than expression of autonomous transposons. To test this, we compiled public RNA-seq datasets from aged human fibroblasts, replicative and drug-induced senescence in human cells, and RNA-seq from aging mice and senescent mouse cells. Indeed, our reanalysis revealed a correlation between transposons expression, intron retention, and transcriptional readthrough across samples and within samples. Both intron retention and readthrough increased with aging or cellular senescence and these transcriptional defects were more pronounced in human samples as compared to those of mice. In support of a causal connection between readthrough and transposon expression, analysis of models showing induced transcriptional readthrough confirmed that they also show elevated transposon expression. Taken together, our data suggest that elevated transposon reads during aging seen in various RNA-seq dataset are concomitant with multiple transcriptional defects. Intron retention and transcriptional readthrough are the most likely explanation for the expression of transposable elements that lack a functional promoter.

Characterization of the novel HLA-DQB1*02:01:01:21Q allele by sequencing-based typing.

HLA-DQB1*02:01:01:21Q differs from HLA-DQB1*02:01:01:01 by one nucleotide substitution in the splice site in the beginning of intron 3.

Developmental and Housekeeping Genes: Two Types of Genetic Organization in the Drosophila Genome.

We developed a procedure for locating genes on Drosophila melanogaster polytene chromosomes and described three types of chromosome structures (gray bands, black bands, and interbands), which differed markedly in morphological and genetic properties. This was reached through the use of our original methods of molecular and genetic analysis, electron microscopy, and bioinformatics data processing. Analysis of the genome-wide distribution of these properties led us to a bioinformatics model of the Drosophila genome organization, in which the genome was divided into two groups of genes. One was constituted by 65, in which the genome was divided into two groups, 62 genes that are expressed in most cell types during life cycle and perform basic cellular functions (the so-called "housekeeping genes"). The other one was made up of 3162 genes that are expressed only at particular stages of development ("developmental genes"). These two groups of genes are so different that we may state that the genome has two types of genetic organization. Different are the timings of their expression, chromatin packaging levels, the composition of activating and deactivating proteins, the sizes of these genes, the lengths of their introns, the organization of the promoter regions of the genes, the locations of origin recognition complexes (ORCs), and DNA replication timings.

Preclinical Development of Antisense Oligonucleotides to Rescue Aberrant Splicing Caused by an Ultrarare ABCA4 Variant in a Child with Early-Onset Stargardt Disease.

Precision medicine is rapidly gaining recognition in the field of (ultra)rare conditions, where only a few individuals in the world are affected. Clinical trial design for a small number of patients is extremely challenging, and for this reason, the development of N-of-1 strategies is explored to accelerate customized therapy design for rare cases. A strong candidate for this approach is Stargardt disease (STGD1), an autosomal recessive macular degeneration characterized by high genetic and phenotypic heterogeneity. STGD1 is caused by pathogenic variants in ABCA4, and amongst them, several deep-intronic variants alter the pre-mRNA splicing process, generally resulting in the insertion of pseudoexons (PEs) into the final transcript. In this study, we describe a 10-year-old girl harboring the unique deep-intronic ABCA4 variant c.6817-713A>G. Clinically, she presents with typical early-onset STGD1 with a high disease symmetry between her two eyes. Molecularly, we designed antisense oligonucleotides (AONs) to block the produced PE insertion. Splicing rescue was assessed in three different in vitro models: HEK293T cells, fibroblasts, and photoreceptor precursor cells, the last two being derived from the patient. Overall, our research is intended to serve as the basis for a personalized N-of-1 AON-based treatment to stop early vision loss in this patient.

On a kneading theory for gene-splicing.

Two well-known facets in protein synthesis in eukaryotic cells are transcription of DNA to pre-RNA in the nucleus and the translation of messenger-RNA (mRNA) to proteins in the cytoplasm. A critical intermediate step is the removal of segments (introns) containing ∼97% of the nucleic-acid sites in pre-RNA and sequential alignment of the retained segments (exons) to form mRNA through a process referred to as splicing. Alternative forms of splicing enrich the proteome while abnormal splicing can enhance the likelihood of a cell developing cancer or other diseases. Mechanisms for splicing and origins of splicing errors are only partially deciphered. Our goal is to determine if rules on splicing can be inferred from data analytics on nucleic-acid sequences. Toward that end, we represent a nucleic-acid site as a point in a plane defined in terms of the anterior and posterior sub-sequences of the site. The "point-set" representation expands analytical approaches, including the use of statistical tools, to characterize genome sequences. It is found that point-sets for exons and introns are visually different, and that the differences can be quantified using a family of generalized moments. We design a machine-learning algorithm that can recognize individual exons or introns with 91% accuracy. Point-set distributions and generalized moments are found to differ between organisms.

Identification of a novel TSC1 gene variant in a patient with atypical vitiligo-like skin lesions: Unveiling the hidden tuberous sclerosis complex.

BACKGROUND: Tuberous sclerosis complex (TSC), an autosomal-dominant disorder, is characterized by hamartomas affecting multiple organ systems. The underlying etiology of TSC is the pathogenic variations of the TSC1 or TSC2 genes. The phenotype variability of TSC could lead to missed diagnosis; therefore, the latest molecular diagnostic criteria for identifying a heterozygous pathogenic variant in either the TSC1 or TSC2 gene filled this gap. Furthermore, the pathogenicity of numerous variants remains unverified, potentially leading to misinterpretations of their functional consequences. METHODS: In this study, a single patient presenting with atypical vitiligo-like skin lesions suspected to have TSC was enrolled. Targeted next-generation sequencing and Sanger sequencing were employed to identify a pathogenic variant. Additionally, a minigene splicing assay was conducted to assess the impact of TSC1 c.1030-2A>T, located in intron 10, on RNA splicing. RESULTS: A novel TSC1: c.1030-2A>T heterozygosis variant was identified in intron 10. In vitro minigene assay revealed that the c.1030-2A>T variant caused exon 11 skipping, resulting in a frameshift in the absence of 112 base pairs of mature messenger RNA and premature termination after 174 base pairs (p.Ala344Glnfs*59). CONCLUSION: The detection of this novel pathogenic TSC1 variant in the patient with atypical vitiligo-like skin lesions enrolled in our study ultimately resulted in the diagnosis of TSC. As a result, our study contributes to expanding the mutational spectrum of the TSC1 gene and refining the genotype-phenotype map of TSC.

Identification of the novel HLA-A*02:01:01:251 allele in a candidate bone marrow donor by next-generation sequencing.

A single nucleotide mismatch within intron 1 differentiates HLA-A*02:01:01:251 from the HLA-A*02:01:01:01 allele.

Comprehensive Evolutionary Analysis of the SMXL Gene Family in Rosaceae: Further Insights into Its Origin, Expansion, Diversification, and Role in Regulating Pear Branching.

SMXL genes constitute a conserved gene family that is ubiquitous in angiosperms and involved in regulating various plant processes, including branching, leaf elongation, and anthocyanin biosynthesis, but little is known about their molecular functions in pear branching. Here, we performed genome-wide identification and investigation of the SMXL genes in 16 angiosperms and analyzed their phylogenetics, structural features, conserved motifs, and expression patterns. In total, 121 SMXLs genes were identified and were classified into four groups. The number of non-redundant SMXL genes in each species varied from 3 (Amborella trichopoda Baill.) to 18 (Glycine max Merr.) and revealed clear gene expansion events over evolutionary history. All the SMXL genes showed conserved structures, containing no more than two introns. Three-dimensional protein structure prediction revealed distinct structures between but similar structures within groups. A quantitative real-time PCR analysis revealed different expressions of 10 SMXL genes from pear branching induced by fruit-thinning treatment. Overall, our study provides a comprehensive investigation of SMXL genes in the Rosaceae family, especially pear. The results offer a reference for understanding the evolutionary history of SMXL genes and provide excellent candidates for studying fruit tree branching regulation, and in facilitating pear pruning and planting strategies.

Distinct functions for the paralogous RBM41 and U11/U12-65K proteins in the minor spliceosome.

Here, we identify RBM41 as a novel unique protein component of the minor spliceosome. RBM41 has no previously recognized cellular function but has been identified as a paralog of U11/U12-65K, a known unique component of the U11/U12 di-snRNP. Both proteins use their highly similar C-terminal RRMs to bind to 3'-terminal stem-loops in U12 and U6atac snRNAs with comparable affinity. Our BioID data indicate that the unique N-terminal domain of RBM41 is necessary for its association with complexes containing DHX8, an RNA helicase, which in the major spliceosome drives the release of mature mRNA from the spliceosome. Consistently, we show that RBM41 associates with excised U12-type intron lariats, is present in the U12 mono-snRNP, and is enriched in Cajal bodies, together suggesting that RBM41 functions in the post-splicing steps of the minor spliceosome assembly/disassembly cycle. This contrasts with U11/U12-65K, which uses its N-terminal region to interact with U11 snRNP during intron recognition. Finally, while RBM41 knockout cells are viable, they show alterations in U12-type 3' splice site usage. Together, our results highlight the role of the 3'-terminal stem-loop of U12 snRNA as a dynamic binding platform for the U11/U12-65K and RBM41 proteins, which function at distinct stages of the assembly/disassembly cycle.

The splicing regulators RBM5 and RBM10 are subunits of the U2 snRNP engaged with intron branch sites on chromatin.

Understanding the mechanisms of pre-mRNA splicing is limited by the technical challenges to examining spliceosomes in vivo. Here, we report the isolation of RNP complexes derived from precatalytic A or B-like spliceosomes solubilized from the chromatin pellet of mammalian cell nuclei. We found that these complexes contain U2 snRNP proteins and a portion of the U2 snRNA bound with protected RNA fragments that precisely map to intronic branch sites across the transcriptome. These U2 complexes also contained the splicing regulators RBM5 and RBM10. We found RBM5 and RBM10 bound to nearly all branch site complexes and not simply those at regulated exons. The deletion of a conserved RBM5/RBM10 peptide sequence, including a zinc finger motif, disrupted U2 interaction and rendered the proteins inactive for the repression of many alternative exons. We propose a model where RBM5 and RBM10 regulate splicing as components of the U2 snRNP complex following branch site base pairing.

Understanding and Rescuing the Splicing Defect Caused by the Frequent ABCA4 Variant c.4253+43G>A Underlying Stargardt Disease.

Pathogenic variants in ABCA4 are the underlying molecular cause of Stargardt disease (STGD1), an autosomal recessive macular dystrophy characterized by a progressive loss of central vision. Among intronic ABCA4 variants, c.4253+43G>A is frequently detected in STGD1 cases and is classified as a hypomorphic allele, generally associated with late-onset cases. This variant was previously reported to alter splicing regulatory sequences, but the splicing outcome is not fully understood yet. In this study, we attempted to better understand its effect on splicing and to rescue the aberrant splicing via antisense oligonucleotides (AONs). Wild-type and c.4253+43G>A variant-harboring maxigene vectors revealed additional skipping events, which were not previously detected upon transfection in HEK293T cells. To restore exon inclusion, we designed a set of 27 AONs targeting either splicing silencer motifs or the variant region and screened these in maxigene-transfected HEK293T cells. Candidate AONs able to promote exon inclusion were selected for further testing in patient-derived photoreceptor precursor cells. Surprisingly, no robust splicing modulation was observed in this model system. Overall, this research helped to adequately characterize the splicing alteration caused by the c.4253+43G>A variant, although future development of AON-mediated exon inclusion therapy for ABCA4 is needed.

Highly Reactive Group I Introns Ubiquitous in Pathogenic Fungi.

Systemic fungal infections are a growing public health threat, and yet viable antifungal drug targets are limited as fungi share a similar proteome with humans. However, features of RNA metabolism and the noncoding transcriptomes in fungi are distinctive. For example, fungi harbor highly structured RNA elements that humans lack, such as self-splicing introns within key housekeeping genes in the mitochondria. However, the location and function of these mitochondrial riboregulatory elements has largely eluded characterization. Here we used an RNA-structure-based bioinformatics pipeline to identify the group I introns interrupting key mitochondrial genes in medically relevant fungi, revealing their fixation within a handful of genetic hotspots and their ubiquitous presence across divergent phylogenies of fungi, including all highest priority pathogens such as Candida albicans, Candida auris, Aspergillus fumigatus and Cryptococcus neoformans. We then biochemically characterized two representative introns from C. albicans and C. auris, demonstrating their exceptionally efficient splicing catalysis relative to previously-characterized group I introns. Indeed, the C. albicans mitochondrial intron displays extremely rapid catalytic turnover, even at ambient temperatures and physiological magnesium ion concentrations. Our results unmask a significant new set of players in the RNA metabolism of pathogenic fungi, suggesting a promising new type of antifungal drug target.

Comparative Mitogenomics Analysis Revealed Evolutionary Divergence among Neopestalotiopsis Species Complex (Fungi: Xylariales).

The genus Neopestalotiopsis consists of obligate parasites that cause ring spot, scab, and leaf blight diseases in higher plant species. We assembled the three complete mitogenomes for the guava fruit ring spot pathogen, Neopestalotiopsis cubana. The mitogenomes are circular, with sizes of 38,666 bp, 33,846 bp, and 32,593 bp. The comparative analyses with Pestalotiopsis fici showed that N. cubana differs greatly from it in the length of the mitogenomes and the number of introns. Moreover, they showed significant differences in the gene content and tRNAs. The two genera showed little difference in gene skewness and codon preference for core protein-coding genes (PCGs). We compared gene sequencing in the mitogenomes of the order Xylariales and found large-scale gene rearrangement events, such as gene translocations and the duplication of tRNAs. N. cubana shows a unique evolutionary position in the phylum Ascomycota constructed in phylogenetic analyses. We also found a more concentrated distribution of evolutionary pressures on the PCGs of Neopestalotiopsis in the phylum Ascomycota and that they are under little selective pressure compared to other species and are subjected to purifying selection. This study explores the evolutionary dynamics of the mitogenomes of Neopestalotiopsis and provides important support for genetic and taxonomic studies.